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1. Introduction 

The formulation of 2D quantum gravity in terms of random matrix variables opened 
the possibility to apply powerful nonperturbative techniques as the method of orthogonal 
polynomials fjj and led to the discovery of unexpected integrable structures associated 
with its different scaling regimes. Originally, a structure related to the KdV hierarchy of 
soliton equations was found in the continuum limit of pure gravity 0. Consequently, the 
partition functions of the A series of models of matter coupled to 2D gravity were identified 
as r-functions of higher reductions of the KP integrable hierarchy [|3| . Similar statement 
concerning the D series was made in The reductions of the KP hierarchy reflect the 
flow structure of the theory of interacting gravity and matter at distances much smaller 
than the correlation length of the matter fields. On the other hand, there are interesting 
phenomena as the topology changing interactions and, in general, all processes involving 
world manifolds with negative curvature, which become important at distances much larger 
than the correlation length of the matter fields. At extra large distances the fluctuations 
of the matter fields are not important and the flow structure of the theory becomes the one 
of pure gravity. Therefore, to study the infrared phenomena, we need a complementary 
description in which the matter is considered as a perturbation of pure gravity. The 
simplest realization of such a description is given by the Kazakov's multicritical points 
of the one- matrix model [pj, where "pure" gravity theories defined as specially weighted 
lattices were shown to be equivalent to theories of gravity coupled to (nonunitary) matter 
fields. 

In this letter we propose a systematic construction of the theories with matter fields as 
deformations of the integrable structure of pure gravity. For this purpose we will exploit 
the microscopic realization of the matter degrees of freedom given by the s/(2)-related 
statistical models: the Ising || and the O(n) f7j models, the SOS and RSOS models and 
their ADE and ADE generalizations @-||. Each of these models can be reformulated as a 
theory of one or several random matrices with interaction of the form tr ln(M a (g>l + l®Mfr). 
Such an interaction can be introduced by means of a differential operator of second order 
H acting on the coupling constants; the partition function of the model is obtained by 
acting with the operator e 11 on the partition function of one or several decoupled one- 
matrix integrals. The Virasoro constraints = for each of the one-matrix integrals 
transform to linear differential constraints e H L^e _H = for the interacting theory. These 
constraints are equivalent to the loop equations, which have been already derived by other 
means || and therefore will not be discussed here. The new point is that the Hirota bilinear 
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equations that hold for each of the one-matrix integrals become, upon replacing the vertex 
operators as V±(z) a — > e H V^j_(,2)e _H , bilinear functional equations for the interacting 
theory. In this way the integrable structure associated with pure gravity is deformed but 
not destroyed by the matter fields. 

We start with deriving the Hirota equations for the one-matrix integral using the 
formalism of orthogonal polynomials. Then we show how these equations generalize for 
the matrix integrals describing theories of matter fields coupled to gravity. We restrict 
ourselves to the microscopic description of the theory, but our construction survives without 
substantial changes in the continuum limit where the loop equations and the bilinear 
equations are obtained as deformations of the Virasoro constraints and the Hirota equations 
in the KdV hierarchy. Their perturbative solution is given by the loop-space Feynman rules 
obtained in . 



2. Orthogonal polynomials and Hirota equations in the one-matrix model 

The partition function of the ensemble of all two-dimensional random lattices is given 
by the hermitian N x N matrix integral 

/oo 
rfMexp (tr^t n M n ) (2.1) 

n=0 

or, in terms of the eigenvalues A;, % = 1, N, 

„ N oo 

Z N [t] = JlldXi exp(^t n A?) Y[(K-Xj) 2 . (2.2) 

i=l n=0 i<,j 

For the time being we restrict the integration in A's to a finite interval [Ax,, A#] on the real 
axis, so that the measure 

oo 

dfit(X) = dXexp (j^nA n ) (2-3) 

n=0 

is integrable for any choice of the coupling constants t n . 



It is known |l l] that the partition function (|2.1|) is a r-function of the KP hierarchy 
of soliton equations. The global form of this hierarchy is given by the Hirota's bilinear 
equations |12| (for a review on the theory of the r-functions see, for example, ||13H -) Below 



we will derive the Hirota equations using the formalism of orthogonal polynomials. 



Introduce for each iV (N = 0, 1, 2, ...) the polynomial 

P N ,t{\) = (det(X-M)) N>t 



1 /■ N 

^- / nd^(Ai)(A-A,) n^-^) 



(2.4) 



It is easy to prove JIJ that the polynomials ( j2.4|) are orthogonal with respect to the measure 
dnt(X)- Indeed, 

Z N J rf/it(Ajv+i)-P/v,t(Ajv+i)A^- +1 
/ IT ^(^O^iv+iCAi, Ajv+i)Ajv(Ai, Ajv)A^ +1 

J i=l 

— - / Yl dfi t (Xi)A N+1 (X u X N+1 ) 



(2.5) 
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N+l 



(— ) N+1 s Xg Ajv+i(Ai, A s , Atv +1 ) 



s=l 



The sum in the integrand is the expansion of the determinant 



Ai • 
A 2 • 


■• Af" 1 
xJV-l 




N+l ■ 


\iV-l 
A 7V+1 
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with respect to its last column. It vanishes for k = 1, N — l, which proves the statement. 
For n = N, one finds 



Hence 



/ ^(A)P 7V , t (A)A Ar = Z N+1 /(N + 1). 



dfjLt(X) P N ,t(X)P k , t (X) = 5 Njh 



(2.6) 



2.7) 



(jV+l)^*]' 

The orthogonality relations (|2.7|) can be written in the form of a contour integral 
namely, 
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-fdz det(z-M) , 



I2.t 



where the integration contour C encloses the point z = and the interval [Ax,, Xr]. Indeed, 
the residue of each of the n poles is equal to the left hand side of ( |2.7| ) multiplied by 
Z k [t]/Z k+1 [t}. 



A set of more powerful identities follow from the fact that the polynomial Ppf t(z) is 
orthogonal to any polynomial of degree less than iV and in particular to the polynomials 
Pk,t' [A], A; = 1, 2, N — 1 where t' = {t' n , n = 1, 2, ...} is another set of coupling constants. 
Written in the form of contour integrals, these orthogonality relations state, for N' < N, 



^, e E" 1 <..-.>^ det(z _ M) ^^_L 



M) / 

' I N't' 



0. 



(2.9) 



The Hirota equations for the KP hierarchy are obtained from eq. (|2.9| ) after expressing 
the mean values of (det ±1 ) in terms of the vertex operators 



oo 

V±(z) = exp ( ± tnZ n ^j exp ( =p In 



n=0 



n=l 



It follows from the definition (|2.2j) that 

Y±{z)-Z N [t] = e ± ^=o tnZn (&et{z-M) ±1 )^Z N [t\ 
and eq. (|2.9| ) is therefore equivalent to 



N,t 



j> dz (v+(z) ■ Z N [t]j (v_(z) ■ Z N >[t'\) = 



(N' < N), 



which is one of the forms of the Hirota equation for KP [|T3| . 
After a change of variables 



f + f 

■En ^ j Vr 



t -t' 



(2.10) 



(2.11) 



(2.12) 



(2.13) 



the Hirota equations (|2.12| ) take its canonical form 



Res 2=0 z N ~ N ' e^n=i 2y ™ zn e £)IT=i ^ ^Z N [x + y]Z N >[x - y] =0 



(2.14) 



where N' < N. The differential equations of the KP hierarchy are obtained by expanding 
( |2. 14j ) in y n . For example, for N' = N, the coefficient in front of yf is 



y=0 



and one finds for the "free energy" u[t] = 24ts log-Z^r the KP equation 



d 2 u d 



du du d 3 u 
~ 4 ^ + 6 ^ + a^J 



0. 



(2.15) 



(2.16) 



The lowest equation in the case N' = N — 1 (modified KP) is the so called Miura trans- 
formation relating the functions u[t] and v[t] = log(Zjv+i[i]/Z/v[i]): 

u = d 2 v - d\v - (d lV ) 2 . (2.17) 



3. Bilinear functional equations in the O(n) model 

The partition function of the O(n) matrix mo del0 is defined by the N x N matrix 
integral J7j 

dMe W \Yt n trM n + - Y 1 , 7- trM n trikH. (3.1) 

L / — ' 2 n + m n\ ml 

n=0 n+m>l 

(the sum over n and m runs the set of nonnegative integers) and describes the ensemble of 
nonintersecting loops on a random graph. A loop of length n is weighted by a factor nT _n 
where T is the temperature of the loop gas. For nG [—2,2], the model exhibits critical 
behavior with spectrum of the central charge C = 1 — 6(g — l) 2 / g , g = — arccos(— n/2); the 
critical behavior for |n| > 2 is this of a branched polymer ||14|| . The matrix integral ( |3.1| ) 
can be again interpreted as the partition function of a Coulomb gas with more complicated 
but still pairwise interaction between particles: 

The interval of integration [Ax,, Xr] should be such that \r < T/2. With the last restriction 
the denominator in the integrand never vanishes. The choice of the integration interval 
does not influence the quasiclassical expansion and hence the geometrical interpretation 
in terms of a gas of loops on the random planar graph. The saddle point spectral density 
is automatically supported by an interval on the half line [— oo, ^] 0. The most natural 
choice for the eigenvalue interval is therefore Ax, — > — oo, Xr — > T/2, with the conditions 

dfi(X) I __ dfj,(X) I _ p. 
d\ \-oo dX If u ' 

The O(n) model reduces to the hermitian matrix model in the limit n— > and/or 
T — > oo, and can be considered as a deformation of the latter in the following sense. Let 
us define the differential operator 



2 



Id 2 ^ 



18 2 ^ T -n-m ( n + TO )| Q Q 



T dti ^ n + m nl ml dt n dt m 

" n+m>l 



(3.3) 



1 We use a Roman letter for the parameter n€ [—2, 2] to avoid confusion with the summation 
index n running the set of natural numbers. 
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acting on the coupling constants. It is easy to see that the partition function ( |3.2|) is 
obtained from the partition function of the one-matrix model by acting with the operator 



3 H. 



zS {n) [t]=tP-Z N [t]. 



(3.4) 



This simple observation will be of crucial importance for our further consideration. It 
means that the integrable structure of the one-matrix model survives in some form in the 
0{n) model. 

The Hirota equations ( |2.12| ) provide, due to the relation (|3.4|), a set of bilinear equa- 
tions for the partition function of the O(n) model 



dz(Y + (z).Z° N ^[t]) (v_(z).Z°^[t'\ 







(N' < N) 



(3.5) 



where the integration contour C_ encloses the interval [Ax,, Xr] and leaves outside point z 
and the interval [T — A^, T — Ax,], and V±(z) = e H V±(z)e _H are the transformed vertex 
operators whose explicit form is 



V± =e ± ^7=1 exp ( t 
(T -2z)- n/2 exp( ±n 



, 1 d ^ 



z dt 
In 



z 



n=l 

1 d 



_d_ 

dt r , 



(T - z) dt 



n=l 



(T-z) 



n 



d_ 

dtr, 



(3.6) 



After being expanded in t n — t' n , eq. ( ^75|) generates a hierarchy of differential equations, 
each of them involving derivatives with respect to an infinite number of "times" t n . 
The functional relations (|3.5|) are equivalent to the bilinear relations 



/ det(^-M) \ /det(T- z - M) n ' 
\det(T - z - MY I „\ det(z-M) ; ~° 



(3.7) 



(N' < N) 



N,t 



[N',t>] 



where ( )nj denotes the average corresponding to the partition function (|3.2|) . These 
relations can be also proved directly by exchanging the order of the integration in the A's 
and the contour integration in z. 
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4. Bilinear functional equations in the ADE and ADE models 

The ADE and ADE matrix models give a nonperturbative microscopic realization 
of the rational string theories with C < 1. Each one of these models is associated with 
a rank r classical simply laced Lie algebra (that is, of type A r , D r , Eqj^) or its afine 
extension, and represents a system of r coupled random matrices. The matrices M a of size 
N a x N a (a = 1, 2, r) are associated with the nodes of the Dynkin diagram of the simply 
laced Lie algebra, the latter being defined by its adjacency matrix G with elements 

Qab _ f 1 if the two nodes are the extremities of a link < ab > ^ ^ 

1 otherwise. 

The partition function ZQ[t\ depends on r sets of coupling constants t = {t^ \ a = 
l,...,r; n = 1,2,...}. The interaction is of nearest-neighbor type and the measure is a 
product of factors associated with the nodes a and the links < ab > of the Dynkin diagram 

/r / r oo 

a=l \a=ln=l 

2^ ^ m + n mini 

a,b m,n=l 

The target space of the A r model is an open chain of r points and its critical regimes 
describe theories of C = 1 — 6 ^ 9 ~g^ matter coupled to gravity, g = 1 ± + 2m, m G 
The target space of the A r -i model is a circle with r + 1 points and its continuum limit 
describes a compactified gaussian field coupled to gravity. The radius of compactification 
is r in a scale where the self-dual radius is r = 2. In this sense the 0(2) model can be 
referred to as the Aq model of the A series. 

Again, the only nontrivial integration is with respect to the eigenvalues X a i (i = 
1, N a ) of the matrices M a : 

m - nn^- w n - a a - L - (4 - 3) 

a=l»=l 1 l<aft> 1 ^ A ™ Ab 0> 

The domain of integration is assumed to be a compact interval [A^, Xr] with \r <T/2. 

The partition function ( f4.3Q can be obtained by acting on the product of r one-matrix 
partition functions ZN a [t a ],a = 1, ...,r, with the exponent of the second-order differential 
operator 



2 

a,b 



, 8 T -n-m l n + m \\ Q Q 

in 1 — — —r- + v 



dtZ dt b n ^ n + m nl ml dt% dtt, 

u u n+m>l n " l 



(4.4) 



7 



namely, 



(4.5) 



a=l 



Each of the one-matrix partition functions on the right hand side of (|4.5| ) satisfies the 
Hirota equations (|2.12| ) . As a consequence, the left hand side satisfies a set of r functional 
equations associated with the nodes a = 1, ...,r of the Dynkin diagram. Introducing the 
transformed vertex operators 



V£(z) =e u \l(z)e- u 



(4.6) 



we find, for each node a = 1, r, 



jf dz(v a + (z) ■ Z%[t\) (v a _(z) ■ Z%[?\) =0 (N' a < N a ) (4.7) 

where the integration contour C_ in ( fl7l) encloses the interval [Al, A^] but leaves outside 
the interval [T — A^, T — AJ . 

Inserting (|4.4| ) in the definition Q4.6| ) we find the explicit form of the dressed vertex 
operators 



H(T-2z)-^ ab e ± ^Zo t > n e W [ T 



-i_d 
dt 



n =l 



dig ^ 

(T-z)" n 9 
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n 



dt b n 



and write the generalized Hirota equations (|4.7| ) in terms of correlation functions of deter- 
minants: 

e Er=i [ ^-*« ]z " / det(z-M a ) 



dz- 



c - U b (T - 2z) Gab \ U b det(T - z - M b ) 



Qab 



N.t 



U b det(T-z-M b 
det(z - M a ) 



(4.9) 



< N a 



[N',t'] 



Here ( ) ^ p denotes the average in the ensemble described by the partition function ( |4.3|) . 
It is possible to prove eq. (|4.9| ) directly by performing the contour integration. It is 
essential for the proof that the interval [T — A^, T — Al] is outside the contour C. 



5. Concluding remarks 

We proposed an alternative nonperturbative formulation of 2D gravity suited for in- 
vestigating the infrared properties of the theory. We used the fact that in the microscopic 
realization based on the si (2) integrable statistical models, the matter can be introduced 
as deformation of pure gravity. As a consequence, the integrable structure of pure gravity 
is not destroyed and manifests itself in the form of bilinear functional equations of Hirota 
type. The main difference between the loop equations and the bilinear equations is that 
the first relate correlation functions boson-like quantities (traces) while the second relate 
correlation functions of fermion-like quantities (determinants). 

The continuum version of these equations will be considered elsewhere. They have 
the same form as the equations in the microscopic theory, the only difference being that 
the vertex operators in the continuum theory are expanded in the half-integer powers of 
the shifted and rescaled variable z. The corresponding coupling constants are related to 
the coupling constants associated with scaling operators as the coefficients in the Laurent 
expansions of the same function at two different points: near the edge of the eigenvalue 
interval and near infinity. 
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